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Abstract 

A novel kind of shape changing (intensity redistribution) collision with potential application to 
signal amplification is identified in the integrable A^-coupled nonlinear Schrodinger (CNLS) equa- 
tions with mixed signs of focusing and defocusing type nonlinearity coefficients. The corresponding 
soliton solutions for N = 2 case are obtained by using Hirota's bilinearization method. The dis- 
tinguishing feature of the mixed sign CNLS equations is that the soliton solutions can both be 
singular and regular. Although the general soliton solution admits singularities we present para- 
metric conditions for which non-singular soliton propagation can occur. The multisoliton solutions 
and a generalization of the results to multicomponent case with arbitrary N are also presented. 
An appealing feature of soliton collision in the present case is that all the components of a soliton 
can simultaneously enhance their amplitudes, which can lead to new kind of amplification process 
without induced noise. 
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I. INTRODUCTION 



It was suggested a long time ago that solitons could be used to carry data at very 
high bit rate in optical communication systems, because of their ability to overcome the 
dispersion limitation through a balance between the self-phase modulation and dispersion 
effects P]. In fact soliton pulses are known to have many other desirable properties, such 
as their robustness against small changes in the pulse shape or amplitude around the exact 
soliton profile leads to treat such changes only as small perturbations on soliton propagation 
3, 0, 1^ . Strictly speaking, the soliton properties can exit only in an ideal fiber. Indeed, in 
a standard telecommunication fiber, the propagation of light pulsesgives rise to a host of 
perturbing effects which inhibit the desirable properties of solitons 5|. One of the strongly 
perturbing effects that comes inevitably into play is the linear attenuation of light along 
the fiber (which is of the order of 0.2dB/km at carrier wavelength 1.55/im), which does not 
permit to keep a constant balance between the self-phase modulation and the group- velocity 
dispersion j^. Although the fundamental soliton propagation cannot be obtained in standard 
fibers, pulse propagation over relatively long distances (and even transoceanic distances) can 
still be obtained through an appropriate combination of dispersion management and optical 
amplification (now mostly based on Er-doped fiber amplifiers and Raman amplifiers) P,!^,!^- 

All the existing amplification processes involve three major ingredients: The first one is a 
pump wave, which serves as a photon reservoir. The second one is an amplification medium, 
that is, a special material in which the pump wave is mixed with the signal. The third 
ingredient is a physical mechanism that can cause a transfer of photons from the pump to 
the signal. Only three types of physical mechanisms have been exploited so far in optical 
amplifiers, namely the laser process used in laser optical amplifiers (e.g. Er-doped fiber 
amplifiers, semi-conductor optical amplifiers) 9'|, the stimulated Raman scattering (used in 
Raman amplifiers) 0, Is] and parametric wave mixing (used in parametric amplifiers) 0,0]. 
Such optical amplifiers do permit to fully compensate the fiber losses, but the amplification 
process is unavoidably accompanied by an undesirable effect of noise generation which is 
commonly referred to as the "amplified spontaneous emission" (ASE) |lO|,lllll]J. Hence, 
one of the most important characteristic parameters of the optical amplifiers developed so 
far is the so-called "noise figure", which serves as a measure of the amount of noise generated 
during the amplification process The ASE increases with the amplifier gain, and there 
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exists an unavoidable amount of noise, known as the amplifier noise figure limit of 3 dB 
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[iS^ . The ASE is one of the major effects that severely degrades the transmission 
quality of ultra-short light pulses over long distances 0, Q, 13 • To radically resolve the 
problem of ASE limitation in high-speed long-distance transmission systems, it is clear that 
the conceptual approach of optical amplification based on the three ingredients mentioned 
above needs to be partially or totally reformulated. 

In the present work, we examine shape changing (intensity redistribution) collisions of 
vector solitons in mixed coupled nonlinear Schrodinger (CNLS) equations, and report some 
results that suggest the possibility of constructing a novel approach of signal amplification. 
The novelty lies in viewing the collision process of solitons as a fundamental physical mech- 
anism for transferring energy from the pump to the signal. The collision involves two vector 
solitons. One of the two solitons, say 5*1, is chosen, to be the signal, while the other soliton 
{S2) serves as the energy reservoir (pump wave). The major virtue of this type of collision- 
based amplification process is that it does not induce any noise, as it does not make use of 
any external amplification medium. 

On the other hand, the study of physical and mathematical aspects of CNLS equations is of 
considerable current interest as these equations arise in diverse areas of science like nonlinear 
optics, optical communication, bio-physics, Bose-Einstein condensates, and plasma physics 
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19|. The fundamental integrable A^-CNLS system is given by the following set 



of equations 



N 



iqj,z + qj,tt + 2fi i^ailqil"^ ] qj = 0, j 



1,2, 



AT, 



(la) 



j=i 



where qj, j = 1,2, . . . , N, is the complex amplitude of the j-th component, the subscripts 
z and t denote the partial derivatives with respect to normalized distance and retarded 
time, respectively, /i represents the strength of nonlinearity (yU > 0) and the coefficients 
cT;'s define the sign of the nonlinearity. System (fTaj) can be classified into three classes as 
focusing, defocusing and mixed types depending on the signs of the nonlinearity coefficients 



(T;'s. The focusing case arises where al 



admits bright soliton solutions 
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o/'s are equal to 1 and the corresponding system 
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These bright solitons are found to 



undergo fascinating shape changing (intensity redistribution) collisions |2l|, |2j, |2^ (for other 




details see for example Refs. [26|, |27|, |28[ ) and such collision properties are not observed in 
systems with defocusing nonlinearity which arises for all a/ = — 1 in (fTaf). The latter system 



possesses either dark solitons in all the components or dark-bright solitons which undergo 



standard elastic collision 



H 29, Q 



A 



defocusing types are given in Refs. 
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so special analytic solutions for the focusing and 
3^ . The third case arises for mixed signs of cxi's 
(that is, +1 or —1). For convenience, we define a^'s for this mixed case as 

ai = 1 for Z = 1, 2, . . . , ra, 

= -1 for I = n + l,n + 2,...,N. (lb) 

Here onwards we refer to Eq. with the above choice of cr^'s as mixed CNLS equations. 
From a physical point of view, system (Q) with N = 2 corresponds to the modified 



Hubbard model in one dimension 



Similar equation, for = 2, is observed in the 



co,.e.. ot e,ecUo.a.„et. pu.e p.pa.at.n i„ ,eft handed .atenal. Q. The above 
of equations ((H) is found to be completely integrable and the corresponding Lax 



pair was obtained in 33|. In their pioneering works Makhankov et al. |33l l35| have shown 
that Eq. (^, for N = 2, admits particular bright-bright, bright-dark, dark-dark type one 
soliton solutions depending upon the asymptotic behaviour of the complex amplitudes qj, 
j = 1,2. Since then very few works have appeared in the literature to analyse the problem 
further 0, |2^ 0, 13 13 i^^^ ^ detailed review of existing results one can refer to 3)- 



Particularly, in a recent work 38], Kanna et al. have obtained stationary solutions of 



mixed CNLS equations with singularities by following an algebraic approach |22l. l4Cll. l4lj . It 
was observed that despite the points of singularities the solutions behave smoothly in finite 
region of the temporal domain. Then the natural question arises as to whether multisoliton 
solutions exhibiting regular behaviour over the entire space-time regions exist and, if so, 
what is the nature of soliton interactions? 

Being motivated by the above fundamental and intriguing aspects, in the present paper 
we perform a detailed study on the bright soliton collision dynamics arising in the mixed 
CNLS system. In particular, we point out that bright solitons of regular type do exist, 
provided the soliton parameters satisfy certain conditions and that the underlying solitons 
undergo novel shape changing/intensity redistribution collisions. The singular solutions turn 
out to be special cases (with specific parametric choices) of the general soliton solutions. An 
important new feature which we identify in the collision process of regular solitons in the 
mixed CNLS case is that after collision a soliton can gain energy in all its components, while 
the opposite takes place in the other soliton. 



This paper is organized as follows. Section II contains the details of Hirota's bilineariza- 
tion procedure j^l for the CNLS equations to obtain soliton solutions. Though the solutions 
obtained in this paper admit both singular and non-singular behaviours, we call them as 
soliton solutions ascribing to their soliton nature in some specific region. In section III, we 
obtain the one and two soliton solutions. Section IV is devoted to a detailed analysis of 
shape changing (intensity redistribution) collisions exhibited by these soliton solutions. The 
procedure to obtain one and two soliton solutions is extended to multisoliton solutions in 
section V. The results of two component case are generalized in a systematic way to the 
multicomponent case with arbitrary number of components following the lines of Ref . j^l • 
Final section is allotted for conclusion. In Appendix A we present the singular station- 
ary three soliton solution for mixed 3-CNLS equations. The multicomponent multisoliton 
solutions of mixed A^-CNLS equations, for arbitrary A^, is given in Appendix B. 



II. BILINEARIZATION OF MIXED CNLS EQUATIONS 



The set of equations ((T)) has been shown to 



certain types of single soliton solutions 



33 



)e completely integrable |33|,|36|, admitting 



for the N = 2 mentioned in 



the Introduction. Here we are concerned with bright-bright multisoliton solutions whose 
intensity profiles vanish asymptotically and with the nature of soliton interactions. 
Let us apply the bilinearizing transformation 



U) 



9 ■ 1 o 
Qj = — , J = 1,2,... 

to Eq. (1) similar to the focusing case o"/ = 1, / = 1, 2, 
set of bilinear equations. 



AT, 



N 



(2) 



This results in the following 



zD, + D^)g^^\f = 0, J = 1,2,..., AT, 



N 



(3a) 
(3b) 



1=1 



where ai is given by Eq. (fTB]) . * denotes the complex conjugate, g^^^^s are complex functions, 
while f{z,t) is a real function and the Hirota's bilinear operators and Df are defined by 



«j™^(|-|.)"(|-|.)'"<.(..w.'.*' 



{z=z',t=t') 



(3c) 
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FIG. 1: Intensity plots of singular one soliton solution of Eq. for N = 2: (a) for the case 
|a[^^| = la^^^l, (b) for the case \oil^\ < |a[^^|. 

The above set of equations can be solved by introducing the following power series expansions 
for g^^^^s and /: 

gU) = ^g'f) + x'g'i^ + J = 1, 2, iV, (4a) 
/ = I + XV2 + XV4 + ..., (4b) 

where x is the formal expansion parameter. The resulting set of equations, after collecting 
the terms with the same power in x, can be solved recursively to obtain the forms of (^''-'^'s 
and /. 

III. SOLITON SOLUTIONS FOR N=2 CASE 

The mixed system (Q) with N = 2 and (Xi = 1, a2 = — 1 is of special physical interest. 
To start with, we consider this particular case. 

A. One soliton solution 

In order to write down the one soliton solution we restrict the power series dl} to the 
lowest order 

g^'^ = X9^^\ J = l,2, f = l + x'f2. (5) 

Then by solving the resulting set of linear partial differential equations recursively, one can 
write down the explicit one soliton solution as 



(6a) 







1 a? \ 


[92 1 







(2) I 1 + (^Vi+Vl+R 



A;i^,sech ( n,n + ^ ] e'^^' , (6b) 



A, V 2 
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FIG. 2: Intensity plots of regular one soliton solution of Eq. ^ ioi N = 2 case. 



where 



r]i = ki(t + ikiz) = r]iii + ir]u,Aj 

I (i)i2 , I (: 
fi [(Ti\ai '1^ + a2\ai 



a?) 



„(2)|2 



, 0-1 



(1 I (Ti|af-'|2 + a-2|«f^P 



-(T2 = 1- 



T7^, J = 1,2, (6c) 



(6d) 



(A;i + klf 

Note that this one sohton solution is characterized by three arbitrary complex parameters 
a"i \ oi{ \ and ki = k^ + iku, where the suffices R and / represent the real and imaginary 
parts, respectively. The quantities ki^Ai and fcii?v42, give the amplitude of the sohton in 
components qi and g2, respectively, subject to the condition 

1 



(Tl\Ai\ +(721^21 



(6e) 



and the soliton velocity in each component is given by 2kii. The position of the soliton is 
found to be 



R 



1 



2k 



IR 



2k 



-In 



IR 



I I (I)l2 , I (2) 12 



(6f) 



From Eq. ()6b|) . it is clear that singular solutions start occurring when = In this 

case, one can easily observe from Eq. ()6d|) that the quantity becomes 0, and one gets the 
solution 



(7) 



which is unbounded. Such an unbounded solution is depicted in Fig. ??(a) for ki = 1 + 
= = 1, and jj, = 1. 
When < becomes negative (so R becomes complex). In this case, singu- 

larity occurs, whenever 



1 - le^le^^i^ = 0, 



(8a) 



or 



ViR 



(8b) 
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Again a singular solution in this case is plotted in Fig. ??(b) for fci = 1 + z, a^^^ = 0.8, 
af^ = 1, and /i = 1. 

However the bright soliton solution is always regular as long as the condition \a^i^\ > 
\a^i^\ is valid in which case is always real and positive, as the denominator (l + e^^'^'^^'^^) 
in Eq. ()(ia|l is always positive definite (as rjin is real) for this choice. This regular one soliton 
solution is shown in Fig. |21for ki = 1 + i, af^ = 1, af'^ = 0.2, and jj, = 1. 

It is also interesting to note here that the polarization vector evolves in a hyperboloid 
defined by the surface l^ip — |^2p = ^ [33], whereas in the Manakov case it is a sphere 
(that is l^ip + |^2p = -) 2J1 • This allows Eq. ffl to admit a rich variety of singular 
and non-singular solutions and makes significant difference in the collision scenario of bright 
solitons arising in the two systems as we will see in the following sections. 



B. Two soliton solution 

To obtain the two soliton solution the power series expansion (j^ is terminated at the 
higher order terms 

^0) ^ X9? + X'9i'\ J = 1,2, (9a) 
/ = l + x'f2 + x'U (9b) 

Then by solving the resultant linear partial differential equations recursively, we can write 
the explicit form of the solution as 



a 



OOg'*?! _|_ (1('^^qV-2 _|_ gr?i+'?i+f?2+<5ij _|_ g»?i+r/2+»?5+52j 



qj = - ^ , J = 1,2, (10a) 

where 

D = I qVi+^I+Ri _|_ qVi+V2+^0 _|_ g^i+»?2+<55 _|_ g»?2+r;J+-R2 _|_ Qrii+Vi+m+V2+R3 



(10b) 



Various quantities found in Eq. (fTUj) . are defined as below: 

Jii _ (^1 ~ ^2)(a?^/«2i - 02^1^11) _ (^2 - ki){ai^\u - ai^K22) 



e 



{ki + kl){kl + k2) {k2 + kl){ki + kl) 



\ki - k2 



12 



= Ti — , tJ'i? , '7*\\i — rT7i^('^ii'^22 - K12K21), (10c) 
{ki + kl){k2 + kl)\ki + kl\^ 
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and 



I (1) (1)* I (2) (2)=! 



hi = 1,2, 



(lOd) 



{k + k*) 

where ai = 1 and (J2 = —1- This solution is characterized by six arbitrary complex param- 
eters a{ \ a{ \ a"2 \ a2'\ ki, and k2- Note that the form of the above two soliton solution 
remains the same as that of the Manakov case (where cxi = +1, cr2 = +1) |^ 13 , except 
for the crucial difference that in the expressions for the parameters k^- in Eq. ()10d|) ai = +1 
and (72 = —1. 

It can also be easily verified that the singular stationary solution for the N = 2 case given 
by Eq. (17) in Ref. |2a] can be obtained for the specific parametric choice 



0, 



0, ku = k2i = 0, fi = 1, 



(11) 



where ?7io and 7720 are two arbitrary real parameters. For this choice of parameters, Eq. (|lUp 
reduces to the form 



where 



D 



1 + 



1 

D 
1 

D 



er,2 _ 



{km - fc2fi)e^^+^^+^2* 

(fcifi - k2R)e^-^+'^'^+^' 
^kl^{kiR + k2R) 



^klR 



^klR 



{kiR - k2Rfe^'+'^l+^^+^*^ 
lQklj^kl^{kiR + k2Rf 



and Tjj is redefined as 



7]j = kjR{t + ikjRz) + r]jo, j = 1, 2, 



(12a) 
(12b) 

(12c) 

(12d) 



where rjjo^s are arbitrary real parameters. The above equation (|12p can be expressed in 
terms of hyperbolic functions as 
'^kiR j kiR + k2R 



Q2 



D 



sinh k2Rt + r]2o + -In 
kiR - k2R \ 2 



kiR - k 



2R 



D 



'2k2R lkiR + k2R . , f . ^ , 1, 

smh kiRt + rjio + -In 



kiR - h 



2R 



Akli^{kiR + k2R)_ 
kiR - k2R 



2 V^kl^{kiR + k2R)_ 



where 



D = -sinh kiRt + k2Rt + r]iQ + r]2o + In 



klR - k2R 



+ I ^ ^"^^ 1 sinh ( kiRt - k2Rt + rjio - 7720 + In 



IR 



k 



2R 



_'2kiRk2R{kiR + k2R] 

k2R 



<'1R 



(13a) 
(13b) 



(13c) 



FIG. 3: Stationary singular two soliton solution for N = 2 case. 



One can check that Eq. (17) given in Ref. 38[ can be re-expressed in terms of hyperbohc 
functions in a form similar to Eq. (|13|). Figure El represents the stationary singular two 
soliton solution at z = for km = 0.2, = —0.25, a^i^ = = —1, af"* = Og^'' = 0, 

and /i = 1. 

Now from the expression ()10|) it can be observed that the denominator can become zero 
for finite values of z and t leading to singular solutions. However, in the case of the general 
two soliton solution (10), it is possible to make the denominator {D in Eq. (10b)) to be 
non-zero for any value of t and z for suitable choice of kj and ctj'^'s, j, I = 1,2. In order to 
do so we rewrite the denominator D (see Eq. ()10b|) ) as 

D = 2e''i«+^2« {e^^^+^^^/'cosh {thr - 7]2R + {Ri - i?2)/2) 



+e °«cos (r/i/ - ?72/ + ^o/) 
+e'^''\os\i {r]iR + r]2R + R3/2) } , (14a) 

where the suffices R and / denote the real and imaginary parts, respectively. Then the 
solution is regular if the above expression is positive for all values of z and t. For this purpose, 
a definite set of criteria can be identified as follows. As in the case of one soliton solution 
in Sec. IIIA, if we choose the parameters af\i,j = 1,2, such that > = 1,2, 

kiR > and k2R > then 

till > 0, K22 > 0. (14b) 

Correspondingly, from Eqs. pOc|) we note that e^^ > and e^'^ > 0, so that e^^''''^^ > 0. 
Then, e^^^~^^'^''^'^cosh(r]iR — r]2R + {Ri — R2) /2) > 0. There is also the other possibility 
Kii < 0, K22 < 0. But it will not lead to regular solution as in this case e^^ and e^^ become 
negative thereby making Ri and R2 complex. 
The term e^^^^ becomes greater than zero if 

«ii/^22 - |/«i2p > 0. (14c) 

Then for this choice e^^/^cosh (rjiR + r]2R + R3/2) is always greater than zero. 

However, the term cos {rju — ri2i + Sqj) oscillates between —1 and 1. So in order that the 
middle term does not compensate the other two terms at any point in space/time resulting 
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FIG. 4: Shape changing (intensity redistribution) colHsion of two solitons in the mixed CNLS 
system for N = 2 case. 



in D being equal to zero, we should have 

Consequently using the expressions ()10c|) in ()14d|) one may deduce the condition 




1 / ft:ii/t22 \kl-k2\ KuK22 - 1/^12^ {ku] , . 



2V^iR^2/? 2\ki + kl\\j kiRk2R \ki + kl\ 

Note that the conditions ()14b|) and p4c|) are necessary conditions to obtain regular solution 
as their falsity will always result in singular solution. Condition p4e|l is a sufficient one as 
its validity con&ms that the solution is always regular. We are unable to prove whether 
condition ()14ej) is also necessary or not due to the complicated form of the function D as 
a function of the variables t and z given by Eq. ()10b|) or ()14ap . It appears that the latter 
can only be checked numerically for given soliton parameter values. In terms of soliton 
parameters the conditions ()14b|) and ()14c|) read as 

|ai'V-|«f^l' > 0, (15a) 
l«2'V-|«ri' > 0, (15b) 

while (ll4eD becomes 



(|aPp-|af)p)(|a«p-|afp) ^ IQkl^kl^ 



Ittf)^^)* - af^af *|2 {km + k2RY + {ku - k2ir 



(15c) 



Thus the two soliton solution satisfying these conditions represent the interaction of two finite 
amplitude bright solitons with definite velocities and their collision behaviour is analysed in 
the foUwing section. 

For illustrative purpose we consider the case kiR > 0, k2R > 0, /i = 1, a^^ = cosh(6'i)e*'^^, 
a'^^ = cosh(6'2)e*'^^, af^ = sinh(6'i)e*'^^, and a*^^ = sinh(6'2)e*'^^, for some arbitrary 6'i,6'2,0i, 
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and 02- Then, the conditions ()14b|) . ()14cp . and ()14ep become 

_ 1 _ 1 

/til — ) /t22 — ; 

2fclR 2fc2R 

|A;i + /c^l^ - 4A;ii?A;2RCosh^ (^12) > 0, 



1 \ki- k2\\ \ki + ^l"^ - AkiRk2RCOsh (6*12) cosh 



+ - V ' 1 > rrriiA:^ (ig) 

4A;i^A;2i?, 4A;iRfc2K|A;i + A;2l' 1/^1 + ^2*1' 

where 6*12 = 61 — 62- A two sohton colhsion process corresponding to the condition (fT^ is 

shown in Fig. |3]for the parameter choice ki = 1.0 + i, ^2 = 1-1 — i,Oi = 0.8, ^2 = 0.2, (pi = 1 

and 02 = 0.3. This colhsion behaviour is analysed in detail in the following section. 

IV. SHAPE CHANGING (INTENSITY REDISTRIBUTION) COLLISIONS OF 
SOLITONS 

Now it is of interest to understand the collision behaviour, shown in Fig. |3J of the regular 
two soliton solution. Figure H] shows the interaction of two solitons Si and 5*2 which are 
well separated before and after collision, in the qi and q2 components. This figure shows 
that after collision, the first soliton Si in the component qi gets enhanced in its amplitude 
while the soliton 5*2 is suppressed. Interestingly, the same kind of changes are observed in 
the second component q2 as well. This collision scenario is entirely different from the one 
observed in the Manakov system where one soliton gets suppressed in one component and 
is enhanced in the other component with commensurate changes in the other soliton. 

On the other hand, conceptually, the collision scenario shown in Fig. EJmay be viewed as 
an amplification process in which the soliton 5*1 represents a signal (or data carrier) while 
the soliton 5*2 represents an energy reservoir (pump). The main virtue of this amplification 
process is that it does not require any external amplification medium and therefore the 
amplification of 5*1 does not induce any noise. 

The understanding of this fascinating collision process can be facilitated by making an 



asymptotic analysis of the two soliton solution as in the Manakov case |21|, |2i |43|. We 
perform the analysis for the choice ki^, k2R > and ku > /c2/- For any other choice the 
analysis is similar. The study shows that due to collision, the amplitudes of the colliding 
solitons 5*1 and 5*2 change from {Al^kiR, Al'kiji) and {Al~k2R, A^~k2R) to {A]^kiR, A^^^kiR) 
and {A\^ k2R, k2R) , respectively. Here the superscripts in A^'s denote the solitons (num- 
ber(l,2)), the subscripts represent the components (number(l,2)) and '±' signs stand for 
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FIG. 5: Elastic collision of two solitons in the mixed CNLS system for N = 2 case. 



±00'. They are defined as 











(aV^ 
[aV/ 




I A]*\ 

\Ar 1 






(17a) 
(17b) 
(17c) 
(17d) 



All the quantities in the above expressions are given in Eq. (Hm) [2ll|24|42|. The analysis 
reveals the fact that, for the non-singular two soliton solution, the colliding solitons change 
their amplitudes in each component according to the conservation equation 

- i^ri' = i^f p - i^ri' = J = 1,2. (18) 

^^ 

This can be easily verified from the actual expressions given in Eq. ()17j) . 

This condition allows the given soliton to experience the same effect in each component 
during collision, which may find potential applications in some physical situations like noise- 
less amplification of a pulse. It can be easily observed from the conservation relation (fTH)) 
that each component of a given soliton experiences the same kind of energy switching dur- 
ing collision process. The other soliton (say 5*2 ) experiences an opposite kind of energy 
switching due to the conservation law 

/oo 
\qj\^dt = constant, j = 1,2, (19) 
'OO 

as required from Eq. (1). 

The asymptotic analysis also results in the following expression relating the intensities of 
solitons 5*1 and S2 in qi and q2 components before and after interaction (see Eq. (fTTj) ). 

\AY?=\Tl\M\-\\ J,/ = 1,2, (20) 
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where the superscripts l± represent the sohtons designated as 5*1 and 5*2 at z — ±00. The 
transition intensities are defined as 



11 -AiA 



2 



i^r = . ":^'^L.. . i = i,2, (21b) 



|l-Ai(a[^'V«2'^)p' 



Ai = — , A2 = — . (21c) 

In fact, this way of energy (amphtude) redistribution can also be expressed in terms of 
inear fractional transformations (LFTs) as in the CNLS system with focusing nonlinearities 

nnn 

2J, |4J, |45| . For example, one can identify from the asymptotic expressions fll7|) that the 
state of 5*1 after interaction (say = 71+ ) is related to its state before interaction (say 
P12 = -jT^) through the following LET, 

A^ 

Ph2 = -nr = ~znrT~~^7T)^ (22a) 



^1 _ '-"11 Pi, 2 + '-"1 2 



where 



Cff = 4^^4^)*(A;2-A;0 + afaf*(fci + A;;), 
CS = -4'^af*(A;2 + fc;), 

CS = afaf>{h - - a«a«*(A;i + fc^). (22b) 

A similar expression can be obtained for soliton 5*2 also. The analysis of such state trans- 
formations preserving the difference of intensities among the components, during collision, 
in the context of optical computing and their advantage in constructing logic gates is kept 
for future study. 

For the standard elastic collision property ascribed to the scalar solitons to occur here 
we need the magnitudes of the transition intensities to be unity which is possible for the 
specific choice 



(1) (2) 



(1) (2) ' 



(23) 



As an example in Fig. |Slwe present the elastic collision for 61 = 62 = 0.2, 0i = 02 = 0.3 (see 
Eq. (fTB|)). with kjS unaltered, j = 1,2, (Note that this choice satisfies the above condition 
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). For all other values of a-"''''s, the soliton energies get exchanged between the solitons 
in both the components as in Fig. |3J 

The other quantities characterizing this collision process, along with this energy redistri- 
bution, are the amplitude dependent phase shifts and change in relative separation distances. 
Their explicit forms can be obtained as in the case of the Manakov model j^, 24|. Explicit 
expressions for the phase shifts $i and $2 of solitons and S2, respectively, during the 
collision are obtained from the asymptotic analysis as 

$1 = = ^^^-^^^-^^\ (24) 

where Ri, R2, and R3 are defined in Eq. (fTUj). 

Then, the change in relative separation distance between the solitons can be expressed 

as 

M2 = tr2-t|2 = %4^$i, (25) 

klRk2R 

where — position of S2 (at z — » ±00) minus position of 5*1 (at z ±00) . 



V. GENERALIZATION OF THE RESULTS TO MULTISOLITON SOLUTIONS 
AND MULTICOMPONENT CASE 

Having discussed the nature of two soliton collision in the two component case (A^ = 2), 
we now wish to study multisoliton collisions for the = 2 as well as > 2 cases. For this 
purpose, we will consider first the three soliton collision scenario for the N = 2 case and 
then extend the analysis to more general cases. 



A. Multisoliton solutions 



It is straightforward to extend the bilinearization procedure of obtaining one and two 



soliton solutions to multisoliton solutions as was done in Ref. 



2J] for the integrable CNLS 



equations with focusing nonlinearity coefficients. Below, we present the form of the three 
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soliton solution for the mixed CNLS equations ((H) as 



where 



D 

D 

g»72 +r?2 +f?3 +»?| +-r3j 



D 



J = 1,2, (26a) 



£) — ]__!_ g»?i+»yl:+-Ri _|_ g»?2+r;5+-R2 _|_ gr)3+»?|+R3 _j_ ^m+Vi+Sio _j_ gf?l'+'?2+<5fo 

_|_g^l+'?3+^20 _|_ QVt+V3+S20 _|_ g'?2+'?|+'530 _j_ g'?2 +^3+5^0 -)- g*?! +»?2+»?2 +^4 
_|_gm+'?S"+'?3+r;5+_R5 _|_ gr)2+»?2+r;3+»?|+-R6 _j_ ^m+vt+m+Va+-riO _|_ g*?! +»?3 +»?5 +^1*0 
_|_g^2+»?3+m+fy|+T-20 _|_ Qm+ri2+^l+V3+r20 _|_ g'?3+»?|+»?l+»?2+™ _)_ g^3+r)3+r)f+»y2+T3*Q 

_^^vi+v*i+V2+V2+m+v*i+R7 _ (26b) 
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Expressions for various quantities given in Eq. ()26|) have the following forms: 



03j 



,^5 



3% 



ki{t + i/cjz), i = 1,2, 3, 

(fci - A;2)(Qi^^^2i - «2^^^ii) _ (fci - fc3)(aS^^/t3i - ag-^^/tii) 



(27a) 



{ki + kl){kl + k2 
{ki - k2){ai^K22 - a^i'^Ku) 



{k^ + kl){kl + kz 



04j 



{k2 - k^){a^2^K:i2 



a 



{ki + k*){k2 + k* 

{ki - k3){a[^\33 - aj^^Ki^) 
(ks + k*){k* + k,) 

{k2 - k3){ai^nu - a3^n2i) 
(kl + k2)ikl + ks) 

{ki - k2){ai^K23 - <y?i^i3) 
{ki + k*){k2 + k*) 

{k2-k{){k3-k^){k3-k2){kl-kl, 
{kl + ki){kl + k2){kl + k3){kl + ki){kl + k2){k 

\f\n2ll^32 - l^22l^3l) + a2''^(/tl2/t31 
{k2 



{k2 + k*){k* + ks 



K.22) 



ik*3 + k2)ik*3 + k3) 

{kl - k3){a[^^K32 - a3^Ki2) 
(ki + k*){k* + k3) 



\kuK22 - 1^121^21) 



X 



2 + h 

^ ^ ki){k3-ki){k3-k2)ik*-kl) 

h ki){kl + k2){kl + k3){kl + ki){kl + A;2)(A;* + ^^3) 

l\l^33l^21 - l^3lf^23) + a:^\n3ll^l3 " 'tll/«33) + «3''H'^23'«11 " '«13't2l) 

(^2 - fci)(fc3 - fci)(A;3 - A;2)(A;3* - kl) 
h){k*2 + k2){kl + k3){kl + fci)(A;3* + k2){kl + 



l \l^22l^33 - /«23/«32) + "2''^ (/^13/^32 " /^33/^12) + (/«12/t23 

m= 1,2,3, e'^^« = ^^i^, e''^" = ^^i^, e^«° : 



{k2 - kl) 


(^2 


-kl) 


{kl + ki){kl + k2) 


{kl + kl){kl + A;2) 


{k3 - kl) 


(k* 


-kl) 


{kl + ki){kl + k3) 




+ ki){k*^ + k3) 


{k3 - k2) 


(k* 


-kl) 


{kl + k2){kl + k3) 


(k* 


+ k2){k3 + k*) 


{k2 - kl) 


(k* 


-kl) 


{k*i + ki){kl + k2) 


(h* 


+ ki){k*3 + k2) 


{kl - /C2) 


(h* 


-kl) 


{kl + ki){kl + k2) 


(k* 


+ ki){kl + k2) 


{k3 - kl) 


(k* 


-kl) 


{kl + ki){kl + /C3) 


m 


+ ki){k* + k3) 



ki + kl' 

Kl't22 



[1^331^11 

[k22'^33 
[1^111^23 
[/t22«13 
[/«33'«12 



^530 

ki + k*' 

K.12l^2l] , 
'^13't3l] , 
'^23 ^32] , 
1^211^13] , 
fi;i2«23] , 
1^131^32] , 



K-221^13) 

1^23 
k2 + k*3 



(27b) 
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\ki - fc2p 




'\k,- 


-fell 


2 




-A;*)(A;3 + A;3*)|A;i + A;*|^ 




'Ik^ + kll^ 



X [(/tll/t22/«33 - /«ll/«23/«32) + (/^12/^23/^31 " /«12/t2l/«33) 

+ {1^21 1^13 1^32 - /^22/^l3/^3l)] , (27c) 

and 

where ai = 1 and (72 = —1. Here a^^\ and al^\ ki, k2 and ^3, j = 1,2,3, are complex 
parameters. 

The solution (j26p also features singular and non-singular behaviours, as in the case of 
one and two soliton solutions depending upon the values of the soliton parameters. Though 
the denominator D in the solution ()26p is cumbersome, possible non-singular conditions can 
be obtained with some effort. Eq. ()26b|l can be rewritten as 



where 



D = 2e''^-+^^«+''^« {e(^^+^«)/2cosh {r],R - r]2R - V3R + {Ri - Re)/^) 

+e(iJ2+iJ5)/2^Qgh iri2R - VlR - V3R + {R2 - i?5)/2) 

j^e^R-M)/\os\i ir],R - 7],n - V2R + (i?3 - i?4)/2) 
^2e^s^oR+rsoR)/2 (cosh(Xi)cos(ri)cos(Zi) - sinh(Xi)sin(ri)sin(Zi)) 
^2e^S20H+r20R)/2 (cosh(X2)cos(r2)cos(Z2) - sinh(X2)sin(F2)sin(Z2)) 
^2e^SsoR+rroR)/2 (cosh(X3)cos(y3)cos(Z3) - sinh(X3)sin(y3)sin(Z3)) 
+e^^/2cosh {r]iR + r]2R + V3R + ^7/2)} , 



^ , (^lOR — TsOi?) ^ , {^20R — T20r) 
Xl = -V3R + ^ , X2 = -ri2R + , 

^ , {Ssor-Tior) (5lO/ + T"30/) 
^3 = -VlR + 2 ' ^ " ~ + 2 ' 

, {S201 + T201) , (^30/ + no/) 

Y2 = r]u - r]si H , Is = r]2i - rj^i + - 



(28a) 



2 

^ ^ (i^ioj - Tsoj) ^ ^ (1^20/ - T20j) ^ ^ (1^30/ - Tioj) ('28b) 
1- 2 '2- 2 '3- 2 ■ ^ > 

Here the suffices R and / denote the real and imaginary parts, respectively. As in the case 
of two soliton solution here also we find the following conditions need to be satisfied for the 
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FIG. 6: Shape changing (intensity redistribution) collision of three solitons in the mixed CNLS 
system for N = 2 case. 



FIG. 7: Elastic collision of three solitons in the mixed CNLS system for N = 2 case. 

solution to be regular: 

e^' >0, i = 1,2,. ..,7, (29a) 

g(i?l+_R6)/2^ g(R2+-R5)/2^ g{i?3+i?4)/2^ gi?,7/2 ^ ^ |g5iofl+T3ofl^ g<520fl+T-20ii^ g^SOfl+nOB | _ 

(29b) 

Note that, the conditions given in ()29a|l are necessary as the falsity of any of them always 
results in singular solution and the last condition (j29b|l is sufficient to ensure that the given 
solution is regular. In fact these conditions can also be expressed in terms of soliton param- 
eters, but due to their cumbersome nature we do not present them here. The appropriate 
choice of parameters can be made by carefully looking at the explicit forms of c^' , e^^° , and 
e""^o,i = 1, ...,7, and j = 1,2,3. 

Such a non-singular solution representing the shape changing (intensity redistribution) 
collision of three solitons 5*1, 5*2, and 5*3 in the two components qi and q2 is shown in Fig. IHl 
for the parameter choice ki = 1 + i, k2 = 1.2 — 0.5i, = 1 — i, jj, = 1, a^^^ = cosh(^i)e*'^i, 
a^^^ = cosh(e2)e''^S a^^^ = cosh{93)e''''\ af^ = sinh(^i)e*'^2 ^ ^^(2) ^ sinh(^2)e''^^ af'^ = 
sinh(6'3)e^'^^ where 9i = 0.8, 02 = 0.4, 63 = 0.2, 0i = 0.5, and 02 = 1-0. From the figure we 
observe that after collision solitons Si and S2 are enhanced in their intensities while there 
occurs suppression of intensity for soliton S3 in both the components qi and q2- It can be 
verified that before and after collision the conservation relation 

i^rr - i^ri' = i^i^r - i^rr = j = 1,2,3, (so) 

/i 

is satisfied, so that the difference of intensities of the solitons between the components qi 
and q2 is preserved before and after the collision process. The standard elastic collision can 
be regained if a^i^ : a'^^ : a^^^ = a^^ : a'^^ : af\ Fig. [7| illustrates such an elastic collision 
for the choice 61=62 = 63 = 0.4, 0i = 02 = 0.5, with same fc/s , j = 1, 2, 3, as in Fig. El 
In a similar manner the four soliton solution can be deduced from Eq.(A2) given in Ref. 



2J| by redefining ^ij cl-S 111 Eq. ()27d|) with i,j running from 1 to 4. We do not present the 



explicit form of it here because of its cumbersome nature. 
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B. Multicomponent case with N>2 



The next step is to generalize the above results for the N = 2 case to arbitrary N with 
N > 2. To do this we follow the earlier work of two of the authors (T.K. and M.L) 2J] on 
the focusing type CNLS equations with all o"/ = 1, / = 1, 2, A^. This study shows that the 
solutions of mixed CNLS equations with N = 2 case can be generalized to arbitrary case 
just by allowing the number of components to run from 2 to and redefining Kj^'s suitably. 

The procedure can be well understood by considering the example of writing down the 
soliton solutions of Eq. for the case = 3. 



1. One soliton solution 



The one soliton solution of mixed 3-CNLS equations obtained by Hirota's method can be 
written as 



a 



a 



(2) 
1 

(3) 



I _|_ QVl+rit+R' 



(31a) 



where 



rji = ki{t + ikiz) 



(31b) 



in which Kh = 

(Ji = 1, 0-2 = 0-3 



M(ai|a«|2+a2|a(2>|2+a3|af'|2l 



ih + kiY 

and without loss of generality we assume either 

1 or o"! = (72 = 1, o"3 = —1. As in the case of A^ = 2, Sec. Ill A, the 
solution is singular if (Ti|af^p + o"2|Q;f''P + o"3|af''p < 0. Otherwise the solution is regular. 
It can be noticed that for any other combination of cr^'s also the above solution satisfies Eq. 
m, for N = 3. 



2. Two soliton solution 



The two soliton solution for the A^ = 3 case is found to possess the same form of Eq. 
(fTn|l . with j = 1, 2, 3, and Kij is given by 



I (1) (1)* I (2) (2)* , (3) (3)^ 
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(32) 



FIG. 8: Stationary singular three soliton solution for = 3 case. 



FIG. 9: Shape changing (intensity redistribution) collision of two solitons in the mixed CNLS 
system, for = 3 case, exhibiting same kind of shape changes for a given soliton in all the three 
components. 

where cr^'s, I = 1,2, 3, can take the value either +1 or —1. Here also the non-singular solution 
exists for the conditions ()14b|) . ()14c|) . and ()14e|) with the redefinition of Kj/s as in Eq. 



3. Three and multisoliton solutions 



A similar analysis can be done for the multisoliton solutions of the multicomponent case 
with arbitrary N . Particularly the three soliton solution of the mixed 3-CNLS equations 
, Eq. with = 3, can be identified to have the form of three soliton solution for the 
N = 2 case with j running from 1 to 3 (that is, now we have three components qi, q2, and 
gs) and here is redefined as 

f (1) (1)* I (2) (2)* , (3) (3)*\ 

i^ij = _^ , «,J = 1,2,3, (33) 

where cr^'s, I = 1,2, 3, can take the value either +1 or —1 (see also Eq. (10) of Ref. 0|). 

It can also be noticed that the stationary singular solution for = 3 case given in Ref. 
38[ results from the above mentioned three soliton solution for the choice 

af) = -e''^o^ =e''^o, = -e'^so^ o^^ = Q, fc,, = 0,/i = l, i^j, ^,j = l,2,3, (34) 

where r/jo's, j = 1, 2, 3, are real parameters. The resulting limiting form reads in terms of 
hyperbolic functions as given in Appendix A. This singular solution at z = is shown in 
Figure IHl The parameters are chosen as a^^"* = —1, a'^^ = 1, a^^^ = — 1, ap'* = 0,i ^ 
j, i,j = 1,2,3, kiR = 0.8, k2R = 0.5, ksR = 0.4 , and /i = 1. 

This procedure can be generalized further to obtain multisoliton solutions of the multi- 
component case with arbitrary A^. For completeness we present the determinant form of the 
A^-soliton solution of A^-component case in Appendix B, following the lines of Ref. for 
the Manakov case. 
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FIG. 10: Shape changing (intensity redistribution) colhsion of two sohtons in the mixed CNLS 
system, for = 3 case, exhibiting same kind of shape changes for a given sohton in the qi and 
components and an exactly opposite colhsion scenario in the q2 component. 

C. Collision scenario in multicomponent cases 

As we increase the number of components the collision behaviour becomes more interest- 
ing. For example, we consider the collision of two solitons in three component (A^ = 3) mixed 
CNLS system. We study the collision dynamics for the following two possible combinations 
of (t's. For illustration, we present two nontrivial scenarios with two different choices of cxi's. 
Case (i): ai = 1, 0-2 = = —1 

For this case, one possible parametric choice for non-singular solution is given by ki = 1.0+i, 

k2 = 0.9 - i, ai^^ = a^2^ = 1 + i, af ^ = 0.2 + 0.4^, ^ = 0.7 + 0.2i, af^ = 0.1 + 0.3i, 

(s) 

0^2 = 0.4 + O.li, and = 1. We plot the two soliton solution corresponding to this 
parameter choice in Fig. IHl The figure shows that after collision there is an enhancement 
(suppression) of intensities (amplitudes) for a given soliton (say soliton Si {S2)) in all the 
three components. Here also one can verify that the difference of intensities is conserved 
according to the conservation law 

|Afr-|4^r-|Arp = -, / = 1,2. (35) 

Case (ii): ai = a2 = 1 (73 = —1 

Next we consider the above possible choice for a's. The nonsingular intensity plots of solitons 
Si and 5*2 are shown in Fig. ^1 The parameters are chosen as ki = 1.0 + 2,^2 = 0.9 — i, 
) = a^^ = M, af ) = 0.2 + 0.4z, = 1, af^ = = 0.3 + 0.2z 

and /i = 1. This figure shows that after collision the intensity of soliton Si{S2) in the first 
and third components gets enhanced (suppressed) while in the second component Si (6*2) is 
suppressed (enhanced) in its intensity. This is a consequence of the conservation given by 
the relation 

|4-|2 + |4-|2-|4-|2 = |A'+|2 + |4+|2_|Af+|2 = i, / = l,2. (36) 
Thus for the two soliton solution of the A^-component case the shape changing (intensity 
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redistribution) collision occurs according to the relation 

N N 



E^d^rP = E^d^M^ = 7, J = 1,2. (37) 



;=i 1=1 ^ 

However the elastic collision occurs for the choice 

(1) (2) (N) 



(1) (2) ••• (N)' 



(38) 



One can also observe that multisoliton solutions for the case N > 2 also undergo the 
above kind of shape changing (intensity redistribution) collisions but with more possible 
ways of energy exchange. 



VI. CONCLUSION 



In this paper we have obtained the bright soliton type solutions of mixed CNLS Eq. 
by applying Hirota's bilinear method. These solutions admit both singular and non- 
singular behaviours depending upon the choice of the soliton parameters. The condition 
for the existence of non-singular one and two soliton solutions for the N = 2 case are 
identified first. Analysing the corresponding collision behaviour reveals the fact that the 
solitons undergo fascinating shape changing (intensity redistribution) collisions with similar 
changes occurring in both components, which is not possible in the well known Manakov 
system. This shape changing (intensity redistribution) collision occurs with a redistribution 
of intensities among the solitons, spread up in two components, in a particular fashion, 
where the intensity difference of the solitons between the two components is preserved after 
collision, and amphtude dependent phase-shifts as well as change in relative separation 
distances also occur. We have extended this study to obtain multicomponent multisoliton 
solutions. Numerical plottings of the solutions show that similar shape changing (intensity 
redistribution) collision behaviour are also observed for the multicomponent case with N > 2 
as in the case of = 2 but with many possible ways of shape variation. Still it is an open 
question to identify the regions in which system (1) admits bright-dark, dark-bright, dark- 
dark soliton solutions. Our study gives an adequate understanding of collision of bright- 
bright solitons arising in system (Q) for mixed signs of nonlinearities. We believe that this 
kind of study will be of interest in the description of magnetic excitations over an anti- 
ferromagnetic vacuum, electromagnetic pulse propagation in left handed materials and so 
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on. In particular one of the most interesting properties of the bright sohtons that we have 
identified in the present work is that the two components of a sohton can be simultaneously 
amplified during a collision process. Using this property, in principle it becomes possible to 
promote the collision process to the rank of a highly efficient amplification process without 
noise generation, in which the gain can be tuned over a relatively large range through a 
careful choice of pre-collision parameters. However, there still remains a lot of work to be 
done to make the fascinating concept of amplifiers with zero noise figure as practical device 
for optical communication systems. For example, an important and challenging issue will 
be to determine whether such amplification process can survive in the presence of strong 
perturbations or in the presence of propagation instabilities. 
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APPENDIX A: SINGULAR STATIONARY THREE SOLITON SOLUTION FOR 
N=3 CASE 



In this appendix we present the singular stationary three soliton solutions of mixed 3- 
CNLS equations. Considering the three soliton solution given by Eq. but now the k^-'s 
are defined as in Eq. (j33|) . the limiting form for the specific choice of parameters given by 
Eq. (jMj) can be deduced as 



or, / / {kiR+k2R){kiR+k3jt) \ 



(l2 



13 



cosh(Ai 



(fc2fl+fc3fl) 



ik2R-k3R) 



cosh(_Bi 



/ / ikiR+k2R)(k2R+k3R) \ 

'^'^2R\I \ (k2R-k^R)(k-iR-k2R) ) 



D 

cosh(A2) 



jkiR+ksR) 
(kiR-ksR) 



cosh(i?2 



^ik^j^z 



nt. I I (kiR+kzR){k2R+k3R) \ 

^'^■iRM \ (k-iR-k^R)(k-iR-k2R) 



D 

sinh(A3) + 
D 



(fcifl+fc2fl) 
(k2R-kiR) 



sinh(S3 



(Ala) 
(Alb) 
(Ale) 
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where 
D 



B2 
As 

Di 
D2 



cosh(D 



{kiR + k2R){kiR + ksn) 



{k2R - kiR){k'iR - kiR) 

(kiR + k2R) {k2R + ksR) , , ^ 



cosh(D2) 

{k2R + ksRiikiR + km) 



{k2R - kiR){k2R - km) 



{k2R + km)t + ^720 + ^30 + 21^^ 

{k2R - km)t + V20 - V30 + ^In 

{kiR + ksR)t + riiQ + 7/30 + ^In 

1 

{kiR - km)t + Vio - V3Q + 
(kiR + k2R)t + ?7io + V20 + ^In 
{kiR - k2R)t + ?7io - r]2Q + ^In 



{.k2R - km){.km - kiR 
{k2R - kiR){km - kiR){ksR - k2R)^ 



cosh(L)4), (Aid) 



IGk^RkmikiR + k2R){kiR + 
{kiR - k2R){kiR + km)klR 


km){k2R + km)"^ 
'^{km-k2R) 


kmikiR + k2R){kiR - km) _ 
{k2R - kiR){km - kiR) 


^^klRklR{kiR + k2R){kiR + 
{kiR - k2R){k2R + km)km 


km)'^{k2R + km) 

^{k3R-k2R) 


k^RikiR + k2R){k2R - km)_ 
{km - kiR){k2R - kiR) 



_lQkl^kl^{kiR + k2R)'^{kiR + km){k2R + km) 
' {km - kiR){k2R + km)klR 
_kfR{kiR + k3R){k3R - k2R)\ ' 



{kiR + k2R + k3R)t + r]io + 7/20 + Vso 

{kiR - k2R){kiR - k3R){k2R - k^R) 
8kiRk2Rk3R{kiR + k2R){kiR + k3R){k2R + k^R 

2{k2R + k3R)k2Rk3R 



{kiR - k2R - k3R)t + r/io - r/20 - ^?30 + In 
{kiR - k2R + k3R)t + 7710 - 7720 + ^30 + In 
{kiR + k2R - k3R)t + ?7io + 7720 - Vso + In 



kiR{k2R — k^R) 
{kzR - kiR)k2R 
'2kiRk3R{kiR + k^R) 
{k2R — kiR)k3R 



(Ale) 



_'2kiRk2R{kiR + /C2i?) 

Particularly, the stationary solution corresponding to the choice given in Eq. (jMj) can be 
easily checked to be the same as the previously reported form given by Eq. (19) in Ref. 
This clearly shows that the more general soliton solutions presented in this paper admit 
singular solutions as special cases which behave as regular and bounded solutions in specific 
regions. 



APPENDIX B: MULTICOMPONENT MULTISOLITON SOLUTIONS 

To write down the multicomponent multisoliton solutions in a formal way we define the 
following (1 X N) row matrix Cg , (A^ x 1) column matrices ipj , 0, j,s = 1,2,. ..,N, and 
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the (A^ X A^) matrix a: 







a 



(2) 



(N) 



( qVi \ 



oV2 



j,s = 1,2,.. .,N, 



^ (7l ... ^ 

as ... 



v 







(Bla) 



where aj, j = 1,2, ...,N, can take value either +1 or —1. Then the A^-sohton solution of 
A^-CNLS system with mixed signs of nonlinearities can be written as 



where 



^9^ 
D ' 

A I (j) 
-I B 
C, 



s = 1,2,3,.. .,N, 



(Bib) 



D 



A I 
-I B 



(Blc) 



in which s denotes the component. Here / is (A^ x A^) unit matrix and the (A^ x A^) matrices 
A and B are defined as 



Bij Kji 



k; ~\~ kj 



1,3 = 1,2,. ..,N, 



(Bid) 



where rji = ki{t + ikiz), ki is complex, f represents the transpose conjugate. Here we remark 
that though presenting the solutions in determinant form seems to be compact, one has to 
explicitly write down the solutions as we have presented in Sees. II - V, for a complete 
characterization and analysis of the solution. This way of expressing the solutions explicitly 
is also useful to identify the particular parameter choice for which the singular stationary 
A^-soliton solution of A^-component case results from the general solutions. In particular, 
by generalizing the Eqs. (jllj) and (j34|) one can identify that the singular stationary A^- 
soliton solution of the A^-component case results from the above solution (jBljl for the choice 
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= (-l)*e^«, i = 1, 2, and a^p = 0, % = 0, = 1, where i j, i, j = 1, 2, 3, N 
and e^'°'s are arbitrary real parameters. 
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FIGURE CAPTIONS 

Fig. 1: Singular one soliton solution of Eq. (^J for N = 2 case. 
Fig. 2: Regular one soliton solution of Eq. ([Q) for = 2 case. 
Fig. 3: Stationary singular two soliton solution for N = 2 case. 

Fig. 4: Shape changing (intensity redistribution) collision of two solitons in the mixed 
CNLS system for N = 2 case. 

Fig. 5: Elastic collision of two solitons in the mixed CNLS system for N = 2 case. 

Fig. 6: Shape changing (intensity redistribution) collision of three solitons in the mixed 
CNLS system for N = 2 case. 

Fig. 7: Elastic collision of three solitons in the mixed CNLS system for the N = 2 case. 

Fig. 8: Stationary singular three soliton solution for = 3 case. 

Fig. 9: Shape changing (intensity redistribution) collision of two solitons in the mixed 
CNLS system, for = 3 case, exhibiting same kind of shape changes for a given soliton in 
all the three components. 

Fig. 10: Shape changing (intensity redistribution) collision of two solitons in the mixed 
CNLS system, for = 3 case, exhibiting same kind of shape changes for a given soliton in 
the gi and gs components and an exactly opposite collision scenario in the q2 component. 
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